Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

Effect of single nucleotide polymorphisms in cytochrome P450 isoenzyme and N-acetyltransferase 2 genes on the metabolism of artemisinin-based combination therapies in malaria patients from Cambodia and Tanzania.

The pharmacogenetics of antimalarial agents are poorly known, although the application of pharmacogenetics might be critical in optimizing treatment. This population pharmacokinetic-pharmacogenetic study aimed at assessing the effects of single nucleotide polymorphisms (SNPs) in cytochrome P450 isoenzyme genes (CYP, namely, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, and CYP3A5) and the N-acetyltransferase 2 gene (NAT2) on the pharmacokinetics of artemisinin-based combination therapies in 150 Tanzanian patients treated with artemether-lumefantrine, 64 Cambodian patients treated with artesunate-mefloquine, and 61 Cambodian patients treated with dihydroartemisinin-piperaquine. The frequency of SNPs varied with the enzyme and the population. Higher frequencies of mutant alleles were found in Cambodians than Tanzanians for CYP2C9*3, CYP2D6*10 (100C → T), CYP3A5*3, NAT2*6, and NAT2*7. In contrast, higher frequencies of mutant alleles were found in Tanzanians for CYP2D6*17 (1023C → T and 2850C → T), CYP3A4*1B, NAT2*5, and NAT2*14. For 8 SNPs, no significant differences in frequencies were observed. In the genetic-based population pharmacokinetic analyses, none of the SNPs improved model fit. This suggests that pharmacogenetic data need not be included in appropriate first-line treatments with the current artemisinin derivatives and quinolines for uncomplicated malaria in specific populations. However, it cannot be ruled out that our results represent isolated findings, and therefore more studies in different populations, ideally with the same artemisinin-based combination therapies, are needed to evaluate the influence of pharmacogenetic factors on the clearance of antimalarials.

Rapid evolution of cancer/testis genes on the X chromosome.

BACKGROUND: Cancer/testis (CT) genes are normally expressed only in germ cells, but can be activated in the cancer state. This unusual property, together with the finding that many CT proteins elicit an antigenic response in cancer patients, has established a role for this class of genes as targets in immunotherapy regimes. Many families of CT genes have been identified in the human genome, but their biological function for the most part remains unclear. While it has been shown that some CT genes are under diversifying selection, this question has not been addressed before for the class as a whole. RESULTS: To shed more light on this interesting group of genes, we exploited the generation of a draft chimpanzee (Pan troglodytes) genomic sequence to examine CT genes in an organism that is closely related to human, and generated a high-quality, manually curated set of human:chimpanzee CT gene alignments. We find that the chimpanzee genome contains homologues to most of the human CT families, and that the genes are located on the same chromosome and at a similar copy number to those in human. Comparison of putative human:chimpanzee orthologues indicates that CT genes located on chromosome X are diverging faster and are undergoing stronger diversifying selection than those on the autosomes or than a set of control genes on either chromosome X or autosomes. CONCLUSION: Given their high level of diversifying selection, we suggest that CT genes are primarily responsible for the observed rapid evolution of protein-coding genes on the X chromosome.

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes.

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

Molecular characterization and phenotypic expression of mutations in genes for gonadotropins and their receptors in humans.

Gonadotropin hormones undergo important dynamic changes during life. Their rise during puberty stimulates gonadal steroid secretion, triggering the development of secondary sexual characteristics and the acquisition of fertility. The full spectrum of possible mutations and polymorphisms in the human gonadotropins and in their receptor genes has been described in recent years. Patients harboring these mutations display a very wide range of phenotypes affecting all aspects of the reproductive axis. An important insight provided by the careful study of these patients lies in the striking gender differences in the phenotypes associated with a given mutation. As a result, the careful study of these rare patients has allowed us to better define the respective roles of luteinizing hormone and follicle-stimulating hormone in normal human pubertal development and in the achievement of full fertility potential in either males or females. In this work, we describe briefly the known mutations in the genes for both gonadotropins and their receptors, and discuss their genotype/phenotype correlations in light of these important gender differences.

Selected polymorphisms in sex hormone-related genes, circulating sex hormones and risk of endometrial cancer.

BACKGROUND: The role of estrogen and progesterone in the development of endometrial cancer is well documented. Few studies have examined the association of genetic variants in sex hormone-related genes with endometrial cancer risk. METHODS: We conducted a case-control study nested within three cohorts to examine the association of endometrial cancer risk with polymorphisms in hormone-related genes among 391 cases (92% postmenopausal at diagnosis) and 712 individually-matched controls. We also examined the association of these polymorphisms with circulating levels of sex hormones and SHBG in a cross-sectional analysis including 596 healthy postmenopausal women at blood donation (controls from this nested case-control study and from a nested case-control study of breast cancer in one of the three cohorts). RESULTS: Adjusting for endometrial cancer risk factors, the A allele of rs4775936 in CYP19 was significantly associated (OR(per allele)=1.22, 95% CI=1.01-1.47, p(trend)=0.04), while the T allele of rs10046 was marginally associated with increased risk of endometrial cancer (OR(per allele)=1.20, 95% CI=0.99-1.45, p(trend)=0.06). PGR rs1042838 was also marginally associated with risk (OR(per allele)=1.25, 95% CI=0.96-1.61, p(trend)=0.09). No significant association was found for the other polymorphisms, i.e. CYP1B1 rs1800440 and rs1056836, UGT1A1 rs8175347, SHBG rs6259 and ESR1 rs2234693. Rs8175347 was significantly associated with postmenopausal levels of estradiol, free estradiol and estrone and rs6259 with SHBG and estradiol. CONCLUSION: Our findings support an association between genetic variants in CYP19, and possibly PGR, and risk of endometrial cancer.

Genome-wide linkage and association analyses to identify genes influencing adiponectin levels: the GEMS Study.

Adiponectin has a variety of metabolic effects on obesity, insulin sensitivity, and atherosclerosis. To identify genes influencing variation in plasma adiponectin levels, we performed genome-wide linkage and association scans of adiponectin in two cohorts of subjects recruited in the Genetic Epidemiology of Metabolic Syndrome Study. The genome-wide linkage scan was conducted in families of Turkish and southern European (TSE, n = 789) and Northern and Western European (NWE, N = 2,280) origin. A whole genome association (WGA) analysis (500K Affymetrix platform) was carried out in a set of unrelated NWE subjects consisting of approximately 1,000 subjects with dyslipidemia and 1,000 overweight subjects with normal lipids. Peak evidence for linkage occurred at chromosome 8p23 in NWE subjects (lod = 3.10) and at chromosome 3q28 near ADIPOQ, the adiponectin structural gene, in TSE subjects (lod = 1.70). In the WGA analysis, the single-nucleotide polymorphisms (SNPs) most strongly associated with adiponectin were rs3774261 and rs6773957 (P < 10(-7)). These two SNPs were in high linkage disequilibrium (r(2) = 0.98) and located within ADIPOQ. Interestingly, our fourth strongest region of association (P < 2 x 10(-5)) was to an SNP within CDH13, whose protein product is a newly identified receptor for high-molecular-weight species of adiponectin. Through WGA analysis, we confirmed previous studies showing SNPs within ADIPOQ to be strongly associated with variation in adiponectin levels and further observed these to have the strongest effects on adiponectin levels throughout the genome. We additionally identified a second gene (CDH13) possibly influencing variation in adiponectin levels. The impact of these SNPs on health and disease has yet to be determined.

Genes mirror geography within Europe.

Understanding the genetic structure of human populations is of fundamental interest to medical, forensic and anthropological sciences. Advances in high-throughput genotyping technology have markedly improved our understanding of global patterns of human genetic variation and suggest the potential to use large samples to uncover variation among closely spaced populations. Here we characterize genetic variation in a sample of 3,000 European individuals genotyped at over half a million variable DNA sites in the human genome. Despite low average levels of genetic differentiation among Europeans, we find a close correspondence between genetic and geographic distances; indeed, a geographical map of Europe arises naturally as an efficient two-dimensional summary of genetic variation in Europeans. The results emphasize that when mapping the genetic basis of a disease phenotype, spurious associations can arise if genetic structure is not properly accounted for. In addition, the results are relevant to the prospects of genetic ancestry testing; an individual's DNA can be used to infer their geographic origin with surprising accuracy-often to within a few hundred kilometres.

Mapping Bias Overestimates Reference Allele Frequencies at the HLA Genes in the 1000 Genomes Project Phase I Data.

Next-generation sequencing (NGS) technologies have become the standard for data generation in studies of population genomics, as the 1000 Genomes Project (1000G). However, these techniques are known to be problematic when applied to highly polymorphic genomic regions, such as the human leukocyte antigen (HLA) genes. Because accurate genotype calls and allele frequency estimations are crucial to population genomics analyses, it is important to assess the reliability of NGS data. Here, we evaluate the reliability of genotype calls and allele frequency estimates of the single-nucleotide polymorphisms (SNPs) reported by 1000G (phase I) at five HLA genes (HLA-A, -B, -C, -DRB1, and -DQB1). We take advantage of the availability of HLA Sanger sequencing of 930 of the 1092 1000G samples and use this as a gold standard to benchmark the 1000G data. We document that 18.6% of SNP genotype calls in HLA genes are incorrect and that allele frequencies are estimated with an error greater than ±0.1 at approximately 25% of the SNPs in HLA genes. We found a bias toward overestimation of reference allele frequency for the 1000G data, indicating mapping bias is an important cause of error in frequency estimation in this dataset. We provide a list of sites that have poor allele frequency estimates and discuss the outcomes of including those sites in different kinds of analyses. Because the HLA region is the most polymorphic in the human genome, our results provide insights into the challenges of using of NGS data at other genomic regions of high diversity.

Sleep in vitro : Erk pathway regulates sleep duration and sleep related genes

To date, for most biological and physiological phenomena, the scientific community has reach a consensus on their related function, except for sleep, which has an undetermined, albeit mystery, function. To further our understanding of sleep function(s), we first focused on the level of complexity at which sleep-like phenomenon can be observed. This lead to the development of an in vitro model. The second approach was to understand the molecular and cellular pathways regulating sleep and wakefulness, using both our in vitro and in vivo models. The third approach (ongoing) is to look across evolution when sleep or wakefulness appears. (1) To address the question as to whether sleep is a cellular property and how this is linked to the entire brain functioning, we developed a model of sleep in vitro by using dissociated primary cortical cultures. We aimed at simulating the major characteristics of sleep and wakefulness in vitro. We have shown that mature cortical cultures display a spontaneous electrical activity similar to sleep. When these cultures are stimulated by waking neurotransmitters, they show a tonic firing activity, similar to wakefulness, but return spontaneously to the "sleep-like" state 24h after stimulation. We have also shown that transcriptional, electrophysiological, and metabolic correlates of sleep and wakefulness can be reliably detected in dissociated cortical cultures. (2) To further understand at which molecular and cellular levels changes between sleep and wakefulness occur, we have used a pharmacological and systematic gene transcription approach in vitro and discovered a major role played by the Erk pathway. Indeed, pharmacological inhibition of this pathway in living animals decreased sleep by 2 hours per day and consolidated both sleep and wakefulness by reducing their fragmentation. (3) Finally, we tried to evaluate the presence of sleep in one of the most primitive species with a neural network. We set up Hydra as a model organism. We hypothesized that sleep as a cellular (neuronal) property may occur with the appearance of the most primitive nervous system. We were able to show that Hydra have periodic rest phases amounting to up to 5 hours per day. In conclusion, our work established an in vitro model to study sleep, discovered one of the major signaling pathways regulating vigilance states, and strongly suggests that sleep is a cellular property highly conserved at the molecular level during evolution. -- Jusqu'à ce jour, la communauté scientifique s'est mise d'accord sur la fonction d'une majorité des processus physiologiques, excepté pour le sommeil. En effet, la fonction du sommeil reste un mystère, et aucun consensus n'est atteint le concernant. Pour mieux comprendre la ou les fonctions du sommeil, (1) nous nous sommes d'abord concentré sur le niveau de complexité auquel un état ressemblant au sommeil peut être observé. Nous avons ainsi développé un modèle du sommeil in vitro, (2) nous avons disséqué les mécanismes moléculaires et cellulaires qui pourraient réguler le sommeil, (3) nous avons cherché à savoir si un état de sommeil peut être trouvé dans l'hydre, l'animal le plus primitif avec un système nerveux. (1) Pour répondre à la question de savoir à quel niveau de complexité apparaît un état de sommeil ou d'éveil, nous avons développé un modèle du sommeil, en utilisant des cellules dissociées de cortex. Nous avons essayé de reproduire les corrélats du sommeil et de l'éveil in vitro. Pour ce faire, nous avons développé des cultures qui montrent les signes électrophysiologiques du sommeil, puis quand stimulées chimiquement passent à un état proche de l'éveil et retournent dans un état de sommeil 24 heures après la stimulation. Notre modèle n'est pas parfait, mais nous avons montré que nous pouvions obtenir les corrélats électrophysiologiques, transcriptionnels et métaboliques du sommeil dans des cellules corticales dissociées. (2) Pour mieux comprendre ce qui se passe au niveau moléculaire et cellulaire durant les différents états de vigilance, nous avons utilisé ce modèle in vitro pour disséquer les différentes voies de signalisation moléculaire. Nous avons donc bloqué pharmacologiquement les voies majeures. Nous avons mis en évidence la voie Erkl/2 qui joue un rôle majeur dans la régulation du sommeil et dans la transcription des gènes qui corrèlent avec le cycle veille-sommeil. En effet, l'inhibition pharmacologique de cette voie chez la souris diminue de 2 heures la quantité du sommeil journalier et consolide l'éveil et le sommeil en diminuant leur fragmentation. (3) Finalement, nous avons cherché la présence du sommeil chez l'Hydre. Pour cela, nous avons étudié le comportement de l'Hydre pendant 24-48h et montrons que des périodes d'inactivité, semblable au sommeil, sont présentes dans cette espèce primitive. L'ensemble de ces travaux indique que le sommeil est une propriété cellulaire, présent chez tout animal avec un système nerveux et régulé par une voie de signalisation phylogénétiquement conservée.